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Three specific t'.sks were outlined in Article I - STATEMENT OF WORK 
of the subject agreement. Of the three, two were accomplished. The 
third was partially accomplished, and will be completed when data from 
a U-2 flight becomes available. Tlie three tasks will be discussed in 
the following paragraphs. 

1.1 A procedure was developed for calculating 24-hour totals of 
evaporation from wet and drying soils that utilizes surface tempera- 
ture data that can be obtained remotely. Its application requires a 
knowledge of the daily solar radiation and the maximum and minimum air 
temperatures (standard Weather Service measurements) . moist surface 
albedo, and maximum and minimum surface temperatures (obtainable from 
surface or airborne sensors) . Details of this procedure are given in 
the appended reprint from SCIENCE 189:991-992, 1975, entitled 
"Estimating evaporation: A technique adaptable to remote sensing," by_ 

S. B. Tdso, R. D. Jackson, and R. J. Reglnato. 

1 ■' Th ' 1 ermal inertia method of remotely sensing soil moisture 
was furi r developed by an experiment in « aich the surface temperatures 
were obtained from thermocouples, hand-held radiation thermometers, and 
the thermal IR band of a multispectral scanner mounted in a NASA air- 
craft. Data wore obtained for both rough and smooth soil surfaces, and 
from soils that v/ere considerably wetter than in previous experiments. 
Results confirm and c.ctend the earlier work on relations between thermal 
inertia and remote sensing, and show that airborne sensors are an 
excellent means of obtaining surface temperatures. Details of this 
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experiment are given in the appended manuscript "Soil water content 
and evaporation determined by thermal parameters obtained from ground- 
based and remote measurements," by R. J. Reginato (USWCL) , S. B. Idso 
(USWCL), J. F. Vedder (NASA/ AMES) , R. T). Jackson (USWCL) , M, B. Blanchard 
(NASA/ AMES), and R. Goettelman (LFE Corp) , which has been accepted for 
publication in the JOURNAL OF GEOPHYSICAL RESEARCH. A second manuscript 
concerning the relationships between reflected solar radiation and soil 
water content is in the rough draft stage. 

A problem that is bothersome to the thermal inertia technique for 
estimating soil moisture is environmental variability. For example, 
a change in the content of water vapor in the atmosphere changes the 
rate of surface cooling at night and this changes the surface tempera- 
ture. A procedure was developed that utilizes Weather Service air 
temperature data to normalize the measured surface temperatures to 
largely account for environmental variability. This procedure is dis- 
cussed in the appended manuscript entitled "Normalization of surface 
temperature data to corapensate for environmental variability in the 
thermal inertia approach to remote sensing of soil moisture," by 
S, B. Idso, R. D. Jackson, and R. J. Reginato. The manuscript has 
been submitted to the JOURNAL OF APPLIED METEOROLOGY. 

1.3 A NASA U-2 aircraft carrying the HCMR simulator was flo^ra 
over Phoenix on 3 September 1975 at 1400 hours. Concurrently, surface 
temperatures were obtained, using a PRT-5 infrared radiation therriiom- 
eter, on 20 ten-acre fields at the Cotton Research Center farm. Bare 
soil temperatures ranged from 59 to 64 C. Alfalfa temperatures were 
33 to 34 C, and cotton plots ranged in temperature from 33 to 37 C. 

Some cotton plots were being irrigated at flight time. This was a 
calibration flight for the HCMR simulator. Processed data from the 
U-2 are not yet available to compare with the ground-based data. 
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Estimating Etaporation: A Technique Adaptable 
to Remote Sensing 

Abslr;u:i. A procciiun' /.s prcsvmeJ lor cakulatina 2-t’lunir luiah ‘>1 evapurauun jmm 
wi itiul clryiii}; xoil.\. Its \tppjic(ition requires a knowledi^e t>l the iluily solar radiation and 
the niaMinwii and inintnumi air temperatures isiandard H eather Service measure- 
menis i, moist surface alhedo f readily estimated or ohiainable Irom a one-lmie measure' 
menu, and maximum and minimum surface temperatures t obtainable from surface or 
airborne sensors). Tests of the technique on a bare field of Avondale h am at Phoenix. 
.‘lri:ona. have shoK/i it to he independent of season. 


liviijiurallort of ualcr from soils uod 
orops is an important I'aoior in managing 
boil) irrigated at/d dryland farming opera- 
tions. It influences the time of seeding, the 
scheduling of irrigations, and various till- 
age practices (/). livaporation is also im- 
portant in determining the water balance 
of watersheds, which allows prediction and 
estimation of runoll' and groundwater re- 
charge, Thus, several techniques have been 
developed over the years to estimate evap- 
oration rates [2). Most of these techniques, 
however, have been of rather limited use- 
fulness in two respects, hirst, they htive de- 
pended on many environmental parame- 
ters and surface characteristics that are 
generally dillicult to measure over extend- 
ed areas, that is, vapor pressure, air tem- 
perature, wind speed gradients, soil wtiter 
content, and surface roughness length. Sec- 
ond, mtiny have been applicable only to 
potential evaporation— the rale that pre- 
vails over a surface of any configuration 
under a given set of meteorological condi- 
tions if there is no saturation deficit at the 
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Fig. 1. Total evaporation induced by thermal ra- 
diation (/.Aj) from a smooth bare Meld of 
Avondale loam at Phoenix, Arizona, as a func- 
tion of the total net thermal radiation (/-■<) eal- 
euiaied as (/Ja /?s) from average valuer, of 7"^ 
and T^. as determined from nighttime data; r ^ 
correlation coefficient. 


surface, that i.s, a eunditiun of nonlimilmg 
water supply (J). 

In light of the seriousness of the current 
and projected world food shortage, wc 
must overcome these limitations and de- 
velop a method of evaporation e.slimaiiog 
readily aJapiahle to rapid application over 
large areas that liandles both the potential 
rale phase of evaporation and the posl- 
poieniiai (falling rate or soil-limiting) 
phase of evaporation, wdiere the surface 
water supply is limiting and aeis to de- 
crease evaporation rates below the poten- 
tial rate that would occur if water were 
nonlimiting. In this repo.n we describe the 
first step in the development of such an 
evaporation estimation technique and its 
initial tests on a field of bare soil, lit addi- 
tion to standard Weather Service measure- 
ments of daily solar radiation and maxi- 
mum and minimum tiir temperatures, it re- 
quires only a one-time r icasuremeni or es- 
timate of moist surface albedo and daily 
measurements of maximum and minimum 
surface temperatures. 

We note first that the evaporation ener- 
gy equivalent {LE) is largely proportional 
to net radiation (/?n) in the potential rate 
phase, and that net radiation can be readily 
subdivided into its two component parts; 
net solar radiation (.yN) and net thermal 
radiation (/-n). Since 5;,< is an external 
forcing function thermally independent of 
evaporation whereas is in part deter- 
mined by the evaporation process by virtue 
of its effects on surface temperature, we as- 
sume that the total 24-hour evaporation is 
directly equal to the daily Sj,.. plus .some 
function of the 24-hour summation of 
that is, wc assume 

LE = Uh+ LEy = Si^+J{Ls) (1) 
where LE^ and LEj arc. respectively, the 
components of the total evaporation in- 
duced by solar and thermal radiation. 

To explicitly derive the relation Z.£‘r = 
we utilized nighttime data, when 
no solar radiation was present. On several 
clear nights we measured evaporation 
from a smooth hare surface of Avondale 


loam witii two weighing Issimcicrs. Wc 
plotted these measurements against night- 
time totals tif /.\, obtained from calcu- 
lations ot (/?A Wi,) where R/^ is the in- 
coming aimosptienc thermal radiation and 
the outgoing surlace thermal radi- 
ation The quantity vvas obtained from 
the Idso-Jackson formula (4) as 

R.s.aff[\ 0.2t,| 

expl 7.77 > 10 T273 7 j-j) (2) 

where n is the Stelaii-liolt/mann constant 
and T.\ is the air temperature measured at 
I m above the surface. The quantity 
was obtained from the Stefan-fiolizmimn 
equation forblackbody radiation as 

Rs^.<if\2 CO 

where 7,, is the surface temperature. Val- 
ues of Iniih 7',\ and were obtained from 
fine-wire, coppcr-consianian thermo- 
couples at 20- or 30-minuie intervals 
through the night. The results (Tig. 1) in- 
dicated that, when the soil surface is moist 
and evaporation is tn the potential rale 
phase. 

LEy = l .'tfi Ls -f 156 t>t) 

In testing our basic hypothesis, wc next 
computed 24-hour reprcseniaiive values of 
Ta and 7s as averages of their maximu.m 
and minimum values and used these 
average values to compute 24-hour totals 
of /-s (which Were all negative). These 
values were then used as the independent 
variable in the linear regression equation 
(Eq. 4) to determine the negative evap- 
oration component to he algebraically 
added to the daily W'e compared 



Fig. 2. Total 24-hour measured evaporation 
(LE) plotted against the 24-hour evaporation 
caieulaied a.s LE Ss -r i 56 + 156. with 

As: in this instance caleuialed from the averages 
of the maximum and minimum values of 7^ and 
7j for the 24-hour period. 
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I-ig. Ratio- of loiiil 24'limir aclual-to-poicntial evaporalion a-, a function of the two thermal 
parameter:; Jetined at the top center o; graph, as detennined for a field of Avondale loam at 
Phoem\, Arizona. 


the results of this procedure with mea- 
surements of potential evaporation (Fig. 
2), which showed that otir basic hypothesis 
produced acceptable results. 

Equations developed to calculated evap- 
oralion will often give good rc.sults in one 
season or climate but not in another (2); 
that is. in a windy, dry situation they may 
do well, but in a calm, humid situation they 
may perform poorly, or vice versa. Wc 
consider our approach potentially adapt- 
able to various situations, for it incorpo- 
rates which is d'~eclly and strongly 
linked to the evaporation rate. For ex- 
ample, between day I and day 3 after 
heavy irrigations of our field (72 by 90 m) 
in February 1962 and again in March 
1971, the average daily wind speed more 
than doubled, greatly increasing the evapo- 
ration rales; yet our calculation procedure, 
which does not explicitly account for wind 
speed, gave equally good results under 
both sets of conditions. Why? Because the 
increased evaporalion rates on the windy 
days lowered the T'j of the soil relative to 
Ta, which resulted in a le.ss negative 7 -n 
flux for the day and a less negative value of 
LEj to be algebraically added to the 
component of the total evaporation. There 
is a similar automatic adjustment for hu- 
midity variations. Over the range of condi- 
tions depicted in Fig. 2, vapor pressure, an- 
other eomponcnl not explicitly accounted 
for in our procedure, varied by a factor of 
5; yet our evaporation calculations were 


equally good over the entire range. In addi- 
tion, no specification of surface type was 
made in developing our technique. Thus, 
we believe that relations similar to the one 
derived for bare soil in Fig. 1 could be de- 
veloped for other surface types such as 
crops. 

For bare soils, wc next confront the 
problem of posipoteniial (falling rale or 
soil-limiting) phase evaporation, where the 
surface becomes dry and evaporation rates 
drop significantly. For this problem we 
again utilized and Idso et al. (J) 
have shown that for several soils, ranging 
from sandy loams to clays, both the maxi- 
mum value minus the minimum value of 
the daily surface soil temperature wave 
( Ts.m< - Ts.min) und the maximum value 
of the surface soil temperature minus the 
air temperature [(7s~'r„)m„ are good 
predictors of soil water pressure potential 
(the w'ork required to move a ^^il mass of 
water against a force field from zero po- 
tential to the point in question), indepen- 
dent of the soil type. Thus, since evapora- 
tion is probably related to water pressure 
potential of the surface soil in the drying 
stages, we fell it would also be related to 
these thermal parameters. 

To lest this idea, we plotted ratios of 
24-hour actual-lo-potenlial evaporation 
against the thermal panimelers K7s„,,„ 
22.5” C] and [(Ts •- 3.5“ C)] 

for several periods after approximate 
lO-cni irrigations of our field (Fig. 3). 


rite potential evaporation for all days was 
taken to be equal to the measured poten- 
tial evaporatiiin at the start of each linie 
series before the surface soil dried, that is. 
tia\ 1 immediatelt after irrigaliun. (SVea- 
iber conditions for all days of each drying 
run were very similar.) On the normali/ed 
basis depicted in Fig. 3, one line adequate- 
ly describes the relation between relative 
evaporation and both of the thernuil pa- 
rameters, Combined with our procedure 
for obtaining actual potential evaporation 
totals for the initial days of such drying 
period!., these results allow estimates of 
actual evaporalion totals to he made 
throughout both the nulential and post- 
potential stages of soil drying, although 
there still remains some uncertainly al the 
transition point between these two re- 
gimes. 

The prime significance of these results 
lies in the fact that they indicate that actual 
evaporalion rates throughout all stages cf 
soil drying may be obtained from remotely 
acquired surface temperatures and routine 
weather network data. Measurements of 
maximum and minimum air temperatures 
are the most basic measurements made at 
all weather stations; solar radiation is rap- 
idly becoming a standard measurement 
also. Moist surface albedo can be obtained 
from information in the literature (6) or 
from a one-time measurement. Thus, max- 
imum and min'mum surface temperatures 
are the only additional data needed for 
successfully estimating evaporalion, and 
these measurements can be made over 
large areas by radiometric means. Such 
temperature measurements may thus be 
capable of specifying actual soil evapora- 
tion rates wherever air temperature and so- 
lar radiation data are available. 

SmfRwooo B. loso 
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Periodic assessment of soil water content and evaporation rates 
can greatly benefit agriculture, hydrology, and civil works [ Idso et al . 
1975a] . Recent research has shown that daily maximum minus minimum 
surface soil temperature and daily maximum surface soil minus air 
temperature (this latter value measured at the time soil maximum occurs) 
can be used to estimate soil-water content and bare soil evaporation 
rates [ Idso et al ., 1975c, d], ihese results were obtained using 
temperatures derived from thermocouples. To extend these techniques I 

j 

to large land areas requires remote radiometric assessment of soil 

10 temperature, whether from just above ground,, from aircraft, or from 

11 satellites. The objective of this paper is to compare measured soil 

12 water contents and bare soil evaporation with estimates derived from 

13 thermal parameters measured by (1) in situ thermocouples, (2) a portable 

14 infrared radiation thermometer just above the ground, and (3) an 

15 aircraft-mounted multispectral scanner. 
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SOIL WATER CONTENT AND EVAPORATION MEASUREMENTS 

2 Three soil moisture conditions were established in a 72 x 90 m 

3 field of bare Avondale loam (fine-loamy, mixed (calcareous), hyperthermic 

Anthropic Torrifluvent) at Phoenix, Arizona. The conditions were: a 

5 continually wet section (#3), a continually dry section (//I), and an 

6 initially wet section (#2) that was allowed to dry during the week of 
the experiment. During the evening of 17 March 1975, two sections 

8 (//2 and #3) were flood-irrigated with about 5 cm of water. A similar 
amount of water was also applied to the surface of a weighing lysiraeter 

10 in section #2. Thereafter, section #3 was replenished with water on the 

11 evenings of 18, 20, and 22 March; section //2 received no additional 

12 water. Section #1 and its associated weighing lysimeter were not I 

13 irrigated at any time during the experiment. 

14 Each section had two surface soil conditions — smooth and rough. 

15 The smooth soil was flat and level and had not been cultivated for the 
10 past 5 years. The lysimeters located within the smooth part of the 

17 section also had smooth surfaces. In contrast, the rough areas were 

18 chiseled and disked to give an uneven surface with roughness elements of 
13 0 to 10 cm. The rough areas were located on the southern portions of all 
20 three sections. 


21 Gravimetric soil moisture samples were taken from all three 

22 moisture treatments on both the smooth and rough plots in identical 

23 fashion to that described previously [ Idso et al . , 1975d] , except 

24 sampling was not continuous but restricted to two 2-hour periods daily: 

25 0430 to 0630 and 1300 to 1500 mountain standard time. Five sets of 
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samples were taken during these time periods (at half-hour intervals) in 
each of the three sections in both the smooth and rough parts. Depth 
intervals sampled were 0 to 0.2, 0 to 0.-5, 0 to 1, 1 to 2, 2 to 4, 

4 to 6, 6 to 8, and 8 to 10 cm. The five samples were averaged to give 
one value per time period for each depth interval for both the smooth 
and rough parts of the three sections. 

Three readings from each of the two lysimeters were obtained 
every 20 minutes, allowing smooth traces of the diurnal soil water 
evaporation trends to be obtained. For comparison with the remote 
sensing techniques, these diurnal trends were Integrated to yielc 24-hou 
totals of bare soil evaporation. Also, daily totals of free water 
evaporation were recorded from buried insulated tanks [ Cooley , 1970], 

To obtain potential soil water evaporation rates for the week's 
experiment (since there was no lysimeter in section 3), the ratio of the 
daily evaporation from the lysimeters during the first 2 days when the 
soil was wet to the daily evaporation from the tanks was computed. This 
ratio (1.04) was then used with tank evaporation data to estimate 
potential soil water evaporation for days when evaporation from the soil 
was below potential. 
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TEMPERATURE MEASUilEMENTS 

Thermocouples » Copper— cons tantan. thermocouples were used to obtair 
air temperature (125 cm above each section) and surface soil temperature 
(about 1 ram below the soil surface) at two locations in the smooth- 
surface parts of the three sections and in the two lysimeters. Tempera- 
tures were recorded every 20 minutes for the week’s experiment. Previous 
research had demonstrated that the minimum surface soil temperature 
occurred just before sunrise, while the maximum occurred about one and 
one-half hours after solar noon. To standardize the measurement times, 
temperatures recorded at 0540 and 0600 hours were averaged for the 
minimum, and readings at 1340 and 1400 hours were averaged for the 
maximum temperature. 

Ground-based radiation thermometer . A portable, precision, 
infrared radiometer (Barnes Engineering^^ PRT-5) was used to measure 
surface soil temperatures at 30-rainute intervals during each of the two 
2-hour sampling periods, from 0430 to 0630 and from 1300 to 1500, for 
the week’s experiment. This instrument yields equivalent blackbody 
temperature with a resolution of + 0.5®C within the 8 to 14 pm band. 
Measured temperatures were corrected for eralttance using a value of 0.96 
for Avondale loam [ Idso and Jackson , 1969). Temperature measurements 
were taken both in the smooth and rough plots of the three sections^, 
and in the two weighing lysimeters. Tlie PRT-5, with a 20® field of 
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^ view, was hand-held at a 1-m height and aimed at a point about 4 ra 
2 distant into the plots, and about 0.4 m high directly over the lysiraeter 
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Airborne scanner . Temperatures were obtained using data from the 

thermal channel (8 to 14 pm wavelength range) of a multispectral scanner 
2 5/ 

(Bendix M S— ) mounted in the NASA NP-3A aircraft. Radiance data in the 

2 

visible region, 0.4 to 1.1 pm, were obtained from the MS. At approxi- 
mately 0540 and 1350 hours each day of the experiment, the aircraft 
passed over the experimental site at about 300 m. The equivalent 
blackbody temperatures were corrected for emittance in the same manner 
as was the PRT-5 data. The effective spatial resolution of the scanner 
for this experiment was about 0.8 x 0.8 ra, with a temperature resolution 
of approximately + 0.5^C. Complete data sets were obtained for 4 of the 
6 days: 18, 19, 20, and 23 March 1975. Instrument malfunction caused 

some data loss on the afternoon of 21 March and the morning of 22 March. 


27 





8 


1 

2 

3 

4 

5 

6 
7 
6 
9 

10 

11 

12 

13 

14 
15| 
16 
17 

I 

18; 

19 

20 
21 
22 

23 

24 

25 

26 


RESULTS AND DISCUSSION 

Previous work [ Idso et al , , 1975b] demonstrated that volumetric 
water content in the upper 2 cm of Avondale loam could be obtained from 
albedo measurements . These measurements appear to be sensitive only to 
the very surface of the soil, whereas surface soil temperatures seem to 
be influenced by soil conditions somewhat deeper. This statement is 
illustrated in Figure 1 which shows a computer-enhanced picture from 
the airborne scanner data of the experimental field on the last day of 
our study. The three sections are as follows: //I on the left was the 

continually dry plot; ffZ in the center was wet initially and allowed to 
dry; and #3 on ib.c right was continually wet. Note that sections 1 and 
2 both appeai^^ light in the visible region, while section 3 is dark. 
However, in the Infrared region there is a definite difference between 
all three sections: #1 being the lightest, if2 intermediate, and #3 the 

darkest. Thus, the thermal infrared region of the spectrum seems to 
hold more promise than does the visible region for the remote assessment 
of soil moisture with depth. 

Soil water content . The relationships between the two thermal 
parameters, daily maximum minus minimum surface soil temperature and 
daily maximum surface soil minus air temperature differential, and 
gravimetric soil water content are shown in Figure 2 for the 0- to ?-cra 
soil depth for both the smooth and rough surface conditions. Earlier 
results had shown a good correlation between the thermal parameters 
and the soil water content in the 0- to 2-cm layer. Data from these 
previous experiments (1970 to 1973) for the 0- to 2-cm depth on smooth 
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Avondale loam for several seasons of the year are shown as solid dots 
in parts A and C. The lines describing the relationships for water 
contents leas than 0,19 were derived from the 1970 to 1973 data using 
linear regression analysis. Also plotted are data obtained in March 
1975 from (1) the in situ surface thermocouples (Q), (2) the ground- 
based Infrared thermometer (Q)* and (3) the airborne multispectral 
scanner (A) , 

In Figure 2A, the recent data for gravimetric water contents less 
than 0.19 are fairly well described by the line derived from the 1970 
to 1973 data. However, at water contents greater than 0.19 the inverse 
relationship does not appear to hold. This is probably explainable as 
follows: for Avondale loam, the soil surface (0- to 2-cm layer) remains 

wet within the water content range from saturation to about 0.19, When 
the surface is wet, the evaporation rate, which greatly influences soil 
temperature, is controlled by meteorological conditions. With uniform 
day-to-day weather, a relatively constant maximum minus minimum soil 
temperature could be expected. For variable meteorological conditions, 

! however, this would not hold true. Thus, the data scatter in Figure 2 

i 

for water contents greater than 0,19 is, at least in part, due to 
variations in meteorological conditions. For water contents below 0.19, 
the rate of moveti»ent of water towards the surface limits evaporation 
and, indirectly, soil temperatures, thereby making evaporation and soil 
temperature less responsive to meteorological conditions. Thus, for 
Avondale loam, less scatter exists below a water content of 0.19, which 
corresponds roughly to the so-called "field capacity," a term used to 
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1 Agreement Is good between the three techniques for measuring 

2 surface soil temperature; In situ thermocouples, ground-based infrared 

3 radiometer, and airborne multispectral scanner. These results are 

4 similar to those obtained by Marlatt {1966]. If water content per se 
5' jure of interest, each soil type or possibly some broader soil 

6 :?.assification unit would be calibrated from rather simple ground-based 

7 aeasurements. Subsequent routine thermal measurements from aircraft or 

8 satellites could be used to characteri?.e the soil moisture status over 

9 Large land areas. 

10 Soil water evaporation . The ratio of actual to potential daily 

11 3oil water evaporation as a function of the two previously described 

12 thermal parameters is shown in Figure 3A and B. Standard error of 

13 sstlmate values were derived in a similar manner as previously described 

1*1 for Figure 2, and are sho;m in Figure 3 as S . The break between the 

y .X 

15 potential and the falling rate stage is at 22“ C for the daily maximum- 
10 mlnlmtom surface soil temperature parameter, and at 3.5“C for the daily 
17 maximum soil minus air temperature parameter. These are the same 
1® values noted in Figure 2 for water content determinations. Intuitively, 
1® one would expect this similarity. When the soil water content is high, 

29 the hydraulic conducting property of the soil is relatively high, thereby 

21 allowing enough water to flow to the soil surface to meet the meteoro- 

22 logical evaporative demand of the potential rate. However, as soil 

23 water becomes limiting the potential evaporation rate cannot be met, 

24 god the relative evaporation rate declines with decreasing water content. 

25 The applicability of these relations to other soils in other areas must 
28 be examined to determine the usefulness of the technique over large 
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areas . 
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Fig. 2. Daily average gravimetric soil water content in the surface 

0- to 2-cm layer of bare Avondale loam daily maximum minus 
minimum surface soil temperature for (A) smooth soil and 
(B) rough soil; and ^ daily maximum soil minus air ' ’ 

temperature for (C) smooth soil and (D) rough soil,- 
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ABSTRACT ' 

A procedure is developed for normalizing surface temperature data ^ 
that are used in the thermal inertia approach to remote sensing of 
soil moisture. The procedure removes data scatter due to environmental 
variability in time and space. Tests of its basic premise on a bare 
soil and a cropped field Indicate it to be conceptually sound. It 
is possible the technique could also be useful in other thermal 
inertia applications, such as lithographic mapping. 
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INTRODUCTION 

A major goal of several scientific groups in the United States 
is to develop a practical procedure for estimating water contents 
near the surfaces of bare soils and throughout the root zones of 
crops from data that can be gathered remotely. Such a feat» if 
accomplished, would open the door to a host of economically Important 
activities, such as predicting world harvests, crop pest outbreaks, 
plant disease epidemics, fertilizer requirements, irrigation needs, 
etc. (Idso, et al. , 1975a). Two basic approaches to achieving this 
goal that have shown substantial Indications of success are to relate 
soil water contents to (1) the magnitudes of the difierences between 
daily maximum and minimum soil or crop canopy temperatures, and 
(2) the differences between maximum soil or crop canopy temperature 
and concurrent air temperature (Idso, et al., 1975b; Idso and 
Ehrler, 1976). 

The first of these proceduxas is what has been known historically 
as the "thermal ifiertia" approach. It has previously been used 
in determining the nature of lunar surface materials prior to 
spacecraft landings (Wesselink, 1948; Jaeger, 1953; Sinton, 1962) and 
in the lithographic mapping of portions of the earth's surface 
(Watson, 1973, 1975; Watson, et al. , 1971; Pohn, et al., 1974; Kahle, 
et al., 1975), based on the fact that the thermal inertia of a 
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given substrate is inversely proportional to the amplitude of its 
diurnal surface temperature oscillation. A problem equally bothersome 
to both of these applications is , environmental variability - the 





non-unlformlty from day-to-day or from season-to-season or from 
place- to-nlace of the external forcing functions of the daily surface 
temperature wave. In this paper we present a solution to this 
problem that may considerably expand the potentials for both remote 
sensing of soil moisture and lithographic mapping. 

THEORY 

The amplitude of the diurnal surface temperature wave of any 
substrate material, be it soil, rock, or plant canopy, is a function 
of both internal and external factors. The internal factors are 
thermal conductivity (X), density (p), and specific heat (C), where 

V = ( 1 ) 

defines what is known as "thermal inertia." The external factors 

include such items as solar radiation, air temperature, atmospheric 

precipitable water extent, cloudiness, wind, aerosol concentration, 

etc. These factors generally are not treated individually in the 

mathematical formalism of thermal inertia analyses, however; but 

their myriad combinations are instead expressed in the single 

resultant forcing function G, the surface heat flux. 

As environmental conditions vary over the earth and in time, G 

may vary considerably, which in turn causes the amplitude of the 

surface temperature wave (0.5 AT ) to vary. This variation is not 

s 

due to variations in P and therefore creates problems for both 





lithographic mapping (based on P discrimination from AT^ measurements) 

2 and soil water content, 0 , estimation (based on 0 vs. AT relations 

V V 8 

3 as surrogate 0^ vs. P relations). 

As a first step in compensating for environmental variability. 


5 we normalize AT^ measurements to what they would have been for some 

6 arbitrary standard value of surface heat flux (G ,). That is, we 

SCu 


8 

9 

10 

11 


transform actual AT data into normalized AT data (AT .. ) via the 

— s s s,Nor 

relationship 



s,Nor. ^std. 


( 2 ) 


12 

13 Thus, in any situation where AT is measured and G is known, we can 

8 

14 transform AT into AT „ , allowing us to make use of a standard 

3 fy f NOIT * 

15 AT ^ vs. P relation that is reasonably Independent of environ- 

s ^ jNoir • 

16 mental conditions. 


17 A problem with this approach is that G is usually not known. 

18 Thus, a surrogate for it must also be found. Air temperature (T^) 

19 would appear to be the ideal candidate for two reasons. First, it is 

20 probably the most commonly measured meteorological parameter on earth. 

21 Second, air temperature respond’ in very similar fashion to the 

22 effects of environmental factors that affect surface temperature. 

23 Indeed, it does so because its diurnal variation is driven by convec- 

24 tive coupling with the surface. Thus, we postulate that 

25 
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and propose that all data be normalized with respect to an 
arbitrary standard diurnal air temperature variation. 


TEST OF THE HEAT FLUX-AIR TEMPERATURE RELATIONSHIP 


During three of our extensive experiments on 0^ vs. AT^ 

relationships in a smooth bare field of Avondale loam (Idso, et al. , 

1975b), we also obtained measurements of soil heat flux at a depth of 

1 cm. These measurements were made with National Instruments 

2 / 

Laboratory— Model HF-1 heat flow discs calibrated by the procedure of 

Idso (1972). Since our analysis of the T^ data indicated that the 

variations in AT^ as 0^ changed were due primarily to changes in 

T and since G . also appeared to be quite invariant, we 

s ) y nxti 

plotted daily vs. AT^ as shoim in Fig. 1, where the data 

were obtained from the nearby National Weather Service Station. The 
results clearly indicate that there is indeed a linear relation 
between ind AT^ of such a nature as to justify equation (3). 


TEST OF THE NORMALIZATION PROCEDURE APPLIED TO BARE SOIL 
Figure 2 contains the original AT^ vs. 0^ data of Idso, et al. 
(1975b), plus some more recent data obtained by Reginato, et al. 

(1976) on the same Avondale loam soil at Phoenix. For each of the 
days represented by data points in Fig. 2, we obtained the maximum 
and minimum air temperatures faeorded by the National Weather Service 


Trade names or company names are included for the benefit of 
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the reader and imply no endorsement or preferential treatment of the 
product listed by the U. S. Department of Agriculture. 



and utilized equation (3) to transform the AT^ values into AT^ 

2 values by arbitrarily assigning aT^ a value of 18 °C. With this 

3 operation the data of Fig. 2 were transformed into the data of Fig. 3, 

4 where the scatter among the data points is seen to be somewhat 

5 reduced .. 

6 The choice of 18 °C foriT^ is completely arbitrary. Any 

7 number could have been chosen. However, to make data from different 

8 locations and seasons compatible, once a number has been chosen, it 

9 must be used exclusively. 

10 

11 TEST OF THE NORMALIZATION PROCEDURE APPLIED TO A CROP 

12 Four separately irrigated 1-hectare plots of Avondale loam planted 

13 to alfalfa at Phoenix, Arizona, were studied from 16 June to 23 July 

14 1975. Every Monday, Wednesday, and Friday, canopy surface tempera- 

15 tures were measured just before sunrise and about an hour and a half 

16 past solar noon. On Tuesdays ar.d Thursdays, only t-e afternoon 

17 measurements were made. The canopy temperatures were measured with a 

2 / 

18 20-degree field-of-viev Barnes PRT-5 infrared thermometer- , hand-held 

19 at about a 45-degree angle with the ground approximately 1 meter 

20 above the crop surface. Preliminary tests using a utility platform 

21 that could be raised 9 meters high indicated that once the alfalfa 
22' was 30 cm high, canopy temperatures did not vary when they were 

23 obtained at viewing angles ranging from 0 to 50 degrees from perpen- 

24 dicular over the height range 1 to 9 meters. 

25 At the same times that canopy temperatures were measured, air 

26 temperatures were measured one meter above the crop canopy by means of 

27 



an aspirated psychroineter . Every Monday, Wednesday, and Friday, we 
also .‘•arapled gravimetric soil water content in each of the four fields 
at 30-cm increments to about 2 meter’s depth. 

The first analyses we made with these data were to test the two 
basic procedures for estimating root-zone soil water contents. Thus, 
in Fig. 4 and 5 we plotted the 1400-hour canopy— air temperature 
differential vs. the volumetric water content of the 0 to 2-meter root 

8 zone, and the 1400-0500-hour canopy temperature differential vs. the 

9 same parameter. Volumetric water contents were obtained by multipli'ing 

10 the measured gravimetric values by the soil's mean bulk density. 

11 The lines drawn on Fig. 4 depict the relation developed by Idso 

12 and Ehrler (1976) for' cotton and sorghum grown on the same soil type. 

13 Our results for alfalfa show essentially the same pattern, where 

14 data for non-water-stressed plants essentially fill up the "bathtub" 

15 part of the graph. The plants we studied were always irrigated at 

16 the proper intervals, however, so that they were never really 

17 stressed. 

18 With this thought in mind, let us consider the data of Fig, 5, 

19 At first glance they appear to be devoid of much meaning. However, it 

20 is noticed that they fall into two major groups: "pre-monsoon" and 

21 "during monsoon." Since our data were all gathered at one location 

22 and we could not traverse great latitude changes to experience 

23 different air temperature regimes due to solar altitude variations, 

24 we conducted our experiment over the period of abrupt climatic change 

25 that occurs with the arrival of Arizona's summer monsoon. During 

26 
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June, Arizona normally experiences very dry weather. However, in 
early July I'c becomes immersed in moist air from both the Gulf of 
Mexico, at high levels, and the Gulf of California, at low levels. 

The low-level source has recently been documented to be the primary 
source (Hales, 1974), which causes the atmospheric precipitable water 
content to about triple in very abrupt fashion. The effect of this 
change in atmospheric humidity is to greatly reduce the amplitude of 
the diurnal air temperature wave, as shown in Fig. 6. Thus, data 
obtained before and after the monsoon's arrival present an ideal 
opportunity for testing our normalization procedure. 

Operating upon the data of Fig. 5, then, in analogous fashion to 
our normalization of the bare soil data that transformed Fig. 2 into 
Fig. 3, we now find Fig. 5 transformed into Fig. 7. The reduction of 
data scatter in this instance is even more than for the bare soil case. 
Indeed, the scatter is cut to only about a third of what it was prior 
to normalization. 

The maximum surface-air temperature differential approach 
cannot claim this same advantage, however, since correct absolute 
values are required for both the surface and air temperatures in order | 
to get a valid differential value. To illustrate, if National 
Weather Service air temperatures are used instead of air temperatures 
measured just above the crop, the plot of Fig. 4 transforms into that 
of Fig. 8. Considerably more scatter is inherent in the data of 
Fig. 8; and the pre-determined soil water content relationship is 
significantly violated. 
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CONCLUDING REMARKS 

In normalizing both the bare soil surface temperature data and 
the alfalfa canopy temperature data, we utilized maximum and minimum 
air temperatures measured at the Phoenix National Weather Service 
Station. Although one cannot expect absolute magnitudes of maximum 
and minimum air temperatures to be the same over a transpiring crop 
or moist soil surface and an asphalt-surrounded airport site several 
kilometers away, the maximum-minimum air temperature differentials 
apparently may be quite similar. This fact greatly increases the 
potential for using the standard thermal inertia approach in remote 
sensing of soil moisture, since no iti situ measurements need to be 
made. It is also possible the technique may be of some usefulness in 
certain lithographic mapping applications, in areas where the 
required air temperatures are available. 
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Fig. 5, 


Fig. 6. 


Fig 


The daily maximum soil heat flux at 1 cm depth in a smooth 
bare field of Avondale loam at Phoenix, Arizona, vs. the 
daily maximum-minimum air temperature differential measured 
at the Phoenix National Weather Service Station. 

The maximum-minimum surface temperature differential of a 
smooth bare field of Avondale loam vs. the average daily 
volumetric soil water content of the uppermost 2 cm. 

Same as Fig, 2, except that the ordinate values of the data 
points have been normalized according to the procedure 
described in the text. 

Maximum canopy-air temperature differential of four different 
fields of mature alfalfa as obtained from measurements made 
at 1400 hours local time vs. the volumetric water content of 
the crops' active root zone. 

The maximum-minimum canopy temperature differential of four 
different fields of mature alfalfa as obtained from measure- 
ments made at 1400 and 0500 hours local time vs. the 
volumetric water content of the crops* active root zone. 

The maximum-minimum air temperature differential obtained 
from official National Weather Service records for Phoenix, 
Arizona, Vi . the mean daily atmospheric precipitable water 
content obtained from National Weather Service dew-point data 
and a procedure outlined by Idso (1969). 

Same as Fig. 5, except that the ordinate values of the data 
points have been normalized according to the procedure 
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described in the text . 





2 


Fig. 8. Same as Fig. A, except that the air temperature data used 
were obtained from the National Weather Service Station, 


3 rather than 1 meter above the crop canopy as in Fig. A. 
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Fig. 1. The daily maximum soil heat flux at 1 cm depth in a smooth bare field of Avondale loam at 
Phoenix, Arizona, vs. the daily maximum-minimum air temperature differential measured at 
the Phoenix National Weather Service Station. 
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Fig. 3. 
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Fig. 4. Maximum canopy-air temperature differential of four different fields of mature alfalfa as obtained 
from measurements made at 1400 hours local time vs. the volumetric water content of the crops' 
active root zone. 
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Fig. 5. The maxitnum-miniTOuni canopy temperature differential of four different fplds of mature alfalfa as 
obtained from measurements made at 1400 and 0500 hours local time vs. the volumetric water content 

of the crops' active root zone. 
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Fig, G. The roaximum-mlniinum air temperature differential obtained from official National Weather Service 
records for Phoenix, Arizona, vs, the mean daily atmospheric precipitable water content obtained 
from National Weather Service dew-point data and a procedure outlined by Idso (1969). 
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Fig. 7. Same as Fig. 5, except that the ordinate values of the data points have been, normalized according 
to the procedure described in the text. 
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